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Abstract 

The inspiral of compact binaries, driven by gravitational-radiation reaction, is investigated 
through 7/2 post-Newtonian (3.5PN) order beyond the quadrupole radiation. We outline the 
derivation of the 3.5PN-accurate binary's center-of-mass energy and emitted gravitational flux. 
The analysis consistently includes the relativistic effects in the binary's equations of motion and 
multipole moments, as well as the contributions of tails, and tails of tails, in the wave zone. How- 
ever the result is not fully determined because of some physical incompleteness, present at the 3PN 
order, of the model of point-particle and the associated Hadamard-type self-field regularization. 
The orbital phase, whose prior knowledge is crucial for searching and analyzing the inspiral signal, 
is computed from the standard energy balance argument. 
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A few years ago it was recognized Q| that improved waveform modelling is crucial to con- 
struct templates for searching and measuring gravitational- waves from inspiralling compact 
binaries with laser- interferometric detectors like LIGO and VIRGO. Since a large number 
of orbital cycles are observable in the frequency band of the detectors, the measurement, 
using the technique of matched-filtering, will be extremely sensitive to those parameters 
that affect the inspiral rate and thus the orbital phase evolution. The orbital phase (which 
is driven by gravitational-radiation reaction) is therefore the crucial quantity to be mon- 
itored for these experiments. Measurement-accuracy analyses have shown that a very 
high post-Newtonian prediction, probably the third post-Newtonian, or even the 3.5PN one 
(i.e. 1/c^), in the case of neutron-star binaries, would be required. Only then the templates 
would be accurate enough over most of the inspiral phase, with reduced cumulative phase 
lags, so that the phasing errors are not significant when one attempts to extract the val- 
ues of the binary's parameters (essentially the masses and spins) from the data. Having in 
hand such high-order post-Newtonian expressions, one could apply resummation methods 
to further improve the convergence of the post-Newtonian series, and make it even more 
accurate for searches as well as parameter estimations Istyi. In this Letter, to provide the 
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essential theoretical input for gravitational- wave data analysis yj, l3 l2l 1^ , we compute the 
orbital phase of compact binaries, both in the time and frequency domains, in the adiabatic 
approximation, at the 3.5PN order. Numerical relativity or approaches such as ^ could de- 
scribe the plunge and merger phases. The latter approach starts from the post-Newtonian 
expansion and goes beyond the adiabatic approximation. Appropriate to the inspiral regime 
Q], we treat the compact bodies as structureless, non-spinning point-particles, moving on 
quasi-circular orbits. Spin effects are known up to 2.5PN order ^ and may be added if 
necessary. 

The first ingredient in the theoretical analysis is the equation of motion of the binary, 
which is used primarily for the calculation of the center-of-mass energy E that is conserved in 
the absence of gravitational-radiation reaction. Recently, the equation of motion of compact 
binaries at the 3PN order has been obtained by means of two different methods, with 
e,uivale,. Ja™* and ScMf. g and Da,.o., Ja.a„ow=M and ScMte. 

employ the Arnowitt-Deser-Misner Hamiltonian formalism of general relativity, while 
Blanchet and Faye , and Andrade, Blanchet and Faye proceed with the post- 

Newtonian iteration of the Einstein field equations in harmonic coordinates. Since the 
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binary's orbit would have been circularized by radiation reaction, the equation of motion is 
of the form 

^ = -^'^^ + ^^-c , (1) 

where = y\ — y2 is the vector separation between the two particles, = dx^ / dt the relative 
velocity, and u the orbital angular frequency {oj = ^tt/P, where P is the period). We denote 
by Fjgg^j. the standard radiation reaction — a resistive force opposite to the relative velocity, 
which arises dominantly at the 2.5PN order. Through 3PN order, the orbital frequency is 
related to the distance r = |x| in harmonic coordinates {via the post-Newtonian parameter 
7 = ^) by a 
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Mass parameters are the total mass m = mi + 1712 and the symmetric mass ratio u = 
77111112/ m? satisfying < z/ < 1/4 (the reduced mass is then = mv). The 3PN coeffi- 
cient depends on two arbitrary constants: a length scale Tq entering the logarithm, and the 
constant A. It was shown in Ref. |p] that rg is merely linked with the choice of harmonic 
coordinates, and has therefore no physical meaning, as it can be eliminated by a change of 
gauge. By contrast, A represents a physical indeterminacy, in the form of a purely numerical 
constant (e.g. a rational fraction), and is probably associated with an incompleteness of the 
Hadamard-type method for regularizing the infinite self-field of point-particles flO], which 
is used to cope with the model of compact objects idealized by Dirac functions (for general, 
non-circular orbits, it is impossible to re-absorb A into a redefinition of the gauge constant 
Tq). The presence of A may be associated with the fact that many integrals composing the 
equation of motion, when taken individually, start depending, from the 3PN order, on the 
internal structure of the bodies, even in the limit where their size tends to zero. However, 
when considering the full equation of motion, we finally expect A to be independent of the 
internal structure of the compact bodies. The constant A is equivalent to the static ambi- 



guity parameter ^static introduced in Refs. P, 0], in the sense that A = — ^cUstatic — l^f^- 
Recently, the value ecstatic = has been obtained by means of a dimensional regularization 
suplementing the ADM-Hamiltonian formalism j^. This result would mean that A = — 
(but we keep A unspecified in this discussion). 
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From now on we shall use in place of the angular frequency uj the dimensionless variable 

X = [^^^S^Y^^ , where / = 2/P = uj/tt is the frequency of the gravitational-wave signal at 

the dominant harmonic. By inverting Eq. Q one finds 7 in terms of the variable x, which 

we shall now consider as an alternative ordering post-Newtonian parameter, 
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As the 3PN equation of motion for general orbits derives from a Lagrangian |ll| (neglecting 
the radiation reaction), one can straightforwardly compute the associated 3PN conserved 
energy. The result, when specialized to circular orbits, reads 
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The good thing to do next is to re-express this energy in terms of the post-Newtonian 
parameter x. Indeed, as x is directly related to the orbital period, the energy will be form 
invariant (the same in different coordinate systems) . We find 0, Q 
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As expected the latter expression is free of the unphysical gauge constant Tq. Since it can be 
checked that for circular orbits there are no terms of order x^*^^, the energy is in fact valid 
up to the 3.5PN order. In the test-mass limit z/ 0, we recover the energy of a particle with 
mass /i in a Schwarzschild background of mass m, i.e. -Etest = l^c^ [(1 — 2x)(l — 3x)^^/^ — l] , 
when developed to the 3.5PN order. 

The second ingredient in this analysis concerns the gravitational wave-form generated 
by the compact binary. More precisely, we need to compute the binary's total energy 
flux at infinity, or gravitational luminosity C, in the post-Newtonian approximation. This 
calculation should take into account the relativistic corrections linked with the description 
of the source (multipole moments), as well as the non-linear effects in the propagation 
of the waves from the source to the far zone. We have applied here a particular wave- 

n n n 

generation formalism |12l . I13L I14| , valid for slowly- moving sources, in which the exterior field 



is parametrized by some specific multipole moments, formally valid to any post-Newtonian 
order Q], and where the observables at infinity are connected to the source moments by 
some non-linear (post-Minkowskian) functional relations, taking into account the various 
effects of tails (see e.g. been speciahzed to the case of 

inspiral wave-forms at the 2.5PN level by Blanchet, Damour and Iyer Furthermore, a 
different formalism, devised by Will and Wiseman Q|, was independently applied to this 
problem and reached equivalent results, reported jointly in Ref. ^|, at the 2PN order. 
The crucial input of any post-Newtonian computation of the flux is the mass quadrupole 
moment (indeed the required post-Newtonian precision of the higher moments is smaller). 
The 3PN quadrupole moment for circular binary orbits, say lij = /i (^Axij + B ^ Vij^ , where 
we neglect a 2.5PN term and denote e.g. ^6ijr^, has been obtained recently by 

Blanchet, Iyer and Joguet fl^, who find the result 
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Note the two types of logarithms entering these formulas at the 3PN order. One type involves 
the same scale Tq as in the equation of motion [see Eqs. (j2))-(jl}]; the other one contains 
a different length scale ro, which is exactly the constant present in the general formalism 
of Refs. 13, [1^ As we know that the constant Tq is pure gauge, it will disappear 
from our physical results at the end. As for tq, it merely represents a convenient scale 
entering the definition of the source multipole moments in Ref. Q], and should cancel out 
when considering the complete multipole expansion of the field exterior to the source. On the 
other hand, besides the harmless constants Tq and ro, there are three unknown dimensionless 
parameters in Eq. ()12|) : ^, k and (. These parameters are analogous to the constant 
A in the equations of motion (see Ref. j3] for their definition in the general case of non- 
circular orbits). They probably reflect an incompleteness of the standard Hadamard self- field 
regularization used in 3]. It is possible that the more sophisticated regularization proposed 



in Ref. Q| could determine some (but maybe not all) of these parameters. However, we 

shall see that, in the case of circular orbits, the energy flux depends only on one combination 

of them: ^ = ^ + 2k + and furthermore that this constant 9 enters the energy flux at 

exactly the same level as A, so that the luminosity given by below depends on one and 

only one combination of 9 and A [to compute the flux one needs the time derivatives of 

the moment ©, and A comes from replacing the accelerations by the equations of motion 

(0)-(j2I)]. More work should be done to determine the values of 9 and A. 

Through 3.5PN order, the result concerning the "instantaneous" part of the total energy 

flux. ,.e. that part which .s generated solely by the multipole moments of the source (not 

counting the tails), is [18|] 
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where 9 = ^ + 2n + C^. The first term represents the Newtonian energy flux coming from the 
usual quadrupole formalism. To the latter instantaneous part of the flux, we must add the 
non- linear tai 



order in Ref. 



effects in the wave zone, which have already been calculated to the 3.5PN 

n 

13| (see Eqs. (5.5a) and (5.9) there). We flnd 
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where C = .577 ■ ■ ■ denotes the Euler constant. What we call here £taii is in fact a com- 
plicated sum of "tails", "tail squares", and "tails of tails", as determined in Ref. Q]. It 
is quite remarquable that so small an effect as a "tail of tail", which constitutes the whole 
3PN coefficient in Eq. (jH)), should be relevant to the present computation, which is aimed 
at preparing the ground for a forthcoming experiment. As we can see, the constant Tq drops 
out from the sum of the instantaneous (jZj) and tail (jH)) contributions — which is normal, 
and constitutes a first test of the calculation. However, the gauge constant Tq does not seem 
to disappear at this stage, but that is simply due to our use in Eqs. dZI)-® of the post- 
Newtonian parameter 7, which depends via the equation of motion on the choice of harmonic 



coordinates. After substituting the frequency-related parameter x in place of 7, i.e. making 
consistent use of the relation Q, we find that Tq does cancel as well — which represents an- 
other test, showing the consistency between the two computations, in harmonic-coordinates, 
of the equation of motion and multipole moments. Finally we obtain 
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The last test (but not the least) is that the expression Q is in perfect agreement, in the 
test-mass limit — > for one of the bodies, with the result following from linear black-hole 
perturbations obtained by Tagoshi and Sasaki 3| (see also Refs. 20|). In particular, the 



rational fraction ^HfUflg^, which is a sum of other fractions appearing in both ((Zj) and 0, 
comes out exactly the same as in the black-hole perturbation theory 19|. 

We shall now deduce the laws of variation of the frequency and phase using a balance 
equation as a fundamental tenet. Namely, we postulate that 



— = -C 
dt 



(10) 



where the binary's gravitational binding energy E is given by Eq. and where the total 
gravitational-radiation luminosity C is the one obtained in Eq. 0. For justifications of the 
validity of the energy balance equation ()1()|1 in post-Newtonian approximations, for either 
point-particle binaries or extended weakly self-gravitating fluids, see Refs. Using the 

previous formulas for E and £, we transform Eq. ()10|) into an ordinary differential equation 
for tu, or, rather, the parameter x. For convenience we adopt a new (dimensionless) time 
variable deflned by 
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where tc denotes the instant of coalescence, at which the frequency tends formally to inflnity 
(evidently, the approximation breaks down well before this point). The solution of the latter 
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differential equation reads 
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Next we compute the binary's instantaneous phase, defined as the angle 0, oriented in the 
sense of the motion, between the separation of the two bodies and, say, the direction of the 
ascending node of the orbit within the plane of the sky. We have ^ = <^ which translates, 
with our notation, into ^ = — ^x^^^, and we can immediately integrate with the result 



/3715 55 
8064 ^ 96^ 



.3/8 _ 



-TTT 



1/4 



+ 



+ 



+ 



+ 



/ 9275495 284875 



+ 



14450688 258048 

831032450749357 



1855 
2048' 



+ 



57682522275840 
123292747421 2255 



53 2 107^ 107, 
— vr^ C H In 

40 56 448 



38645 
172032 

r 



+ 



65 
2048" 



1/ 7rln 



256 



4161798144 
/ 188516689 4J 



2048 



-TT 



385 

Is"' 



^25 



141769 



-u 



-v 



55 
— / 
16 



TTT 



154565 
1835008' 



1179625 
1769472' 



r 



-1/8 



-1/4 



(13) 



173408256 516096 516096 
The constant Tq is related to a constant phase that is simply fixed by the initial conditions 
when the frequency of the wave enters the detector's bandwidth. Finally it can be useful to 
dispose of the expression of the phase in terms of the frequency x. For this we have 
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TABLE I: Contributions to the accumulated number N = :^('^isco "^seismic) of gravitational- wave 
cycles. Frequency entering the bandwidth is /seismic = 10 Hz; terminal frequency is assumed to 

3 

be at the Schwarzschild innermost stable circular orbit /isco = ^i/i^^Qra • "^^^ 3PN term depends 
on the unknown parameter 9 = 9 — (we have 9 = 9 + using the value of A following from 
"^static = 0)- 
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where Xq is determined by the initial conditions (like tq). As a rough estimate of the relative 
importance of each of the various post-Newtonian terms for LIGO/VIRGO-type detectors, 
we give in Table I their contributions to the accumulated number M of gravitational-wave 
cycles (see also Table I in Ref. for the contributions of spin-orbit and spin-spin effects). 
The result for Aspn is given as a function of the combination of parameters 6 = 6 —'^X that 
enters Eq. p4j) . As we can see, if 9 is of the order of unity, we reach with the 3PN or 3.5PN 
approximations an acceptable level of, say, a few cycles, that roughly corresponds to the 
demand which was made by data-analysists in the case of neutron-star binaries P,|^. Indeed, 
the above estimation suggests that the neglected 4PN terms will yield some systematic errors 
that are, at most, of the same order of magnitude, i.e. a few cycles, and perhaps much less. 
However, this conclusion is quite sensitive to the exact value of 9. If 9 is of the order of 
10, we find that the 3PN term is nearly as important numerically as the previous 2PN 
approximation. Finally, in order to define the theoretical template of the compact binary 
inspiral, one should insert the previous 3.5PN-accurate expressions of the frequency and 
phase into the two polarisation wave-forms and h^. A standard practice is to neglect in 
hj^ and /ix all the harmonics but the dominant one / at twice the orbital frequency (i.e. the 
so-called restricted post-Newtonian approximation), but it is better to define the template 
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by means of the 2PN-accurate polarization wave- forms calculated in Ref. j23[. 
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